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Abstract. We show that that vector field-based models of the ether generi- 
cally do not have a Hamiltonian that is bounded from below in a fiat spacetime. 
We also demonstrate that these models possess multiple light cones in flat or 
curved spacetime, and that the non-lightlike characteristic is associated with 
an ether degree of freedom that will tend to form shocks. Since the field 
equations (and propagation speed) of this mode is singular when the timelike 
component of the ether vector field vanishes, we demonstrate that linearized 
analyses about such configurations cannot be trusted to produce robust ap- 
proximations to the theory. 



1. Introduction 

Recently there has been interest in models in which constants of nature are 
elevated to the status of fields |], |j| . This is partly due to recent experimental 
evidence for a variation of the fine structure constant with time [Q, El pj, 0] , and 
partly from theoretical interest in, for example, solving the horizon and flatness 
problems in cosmology |^|, [jj or modeling Planck scale deviations from flat spacetime 
field theory in the very early universe 

Along with these models, either explicitly or implicitly, comes a breakdown of 
Lorentz invariance. This appears as contributions to the action that are explicitly 
written in a preferred frame [jlO|, or in a more subtle manner where, although no 
preferred frame is chosen, one or more of the equivalence principles are violated [ jTl[ 
pjl , |J| . Indeed, the motivation for these models is to embody the idea of a "varying 
speed of light" , which should conflict with Einstein's postulate of the invariance of 
the propagation velocity of light. 

Among the latter class of models we would include those that have similar prior- 
geometric structure jl3| , the so-called "fc-inflation" and "fc-essence" models [|l4| |l5| , 
and vector field models of the ether |l6|, [l7| [l8| . Although these models are not nec- 
essarily intended introduce fields with dynamically determined light velocities, this 
nevertheless results from their construction. In the same way, higher-dimensional 
models Jl|, |o[ |l| also generically end up with multiple light cones as a result of 
dimensional reduction. 

The class of model that we want to discuss here are those considered, for example, 
in Jl(| [It], |l8). These models use a (co-) vector field to model the preferred 
frame representing the ether, introducing a constraint in the action via a Lagrange 
multiplier A that forces to be a timelike, null or spacelike vector. The model is 
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described by the action: 

S=-Jd^R + Jd t i(-^P 2 + ^X(A 2 -l)+A li Ji^, (1) 

where we are using a metric with signature (+, — , — , — ), have chosen 16t:G = 1 = c, 
and have normalized the vector field to remove a possible arbitrary coefficient 
multiplying the F 2 term. We have also written d/i — yj—gd^x, A 2 = g^ v A^A v , 
and F 2 — g^ a g vl3 F^F a p. The source-coupling term A^J 11 is included to represent 
the possible non-derivative coupling of the ether to matter fields through a matter 
action: 5 ma tter(g, A,...), and so one may think of = 9S , ma ttcr/9A M . 

There is a fairly long history of difficulties associated with coupling vector fields 
to gravity. Kuchaf [^2| has noted that in general, fields that are derivatively- 
coupled to gravity may propagate off the light cone as determined by the spacetime 
metric. Examples of this happening when 'higher spin' fields were coupled to gravity 
have been in existence for quite some time p3| , p4| |25| , |26|, and linear models 
coupling a vector field to gravity are known to possess instabilities unless they are 
one of Maxwell, Proca, or a purely longitudinal vector field, minimally-coupled to 
gravity ||27|j . In this paper we demonstrate that the models described by ([[]) suffer 
from similar problems, not so much to point out flaws in the literature as to further 
illustrate the consequences of modifying the structure of vector field models. As 
these issues are serious, they play a complimentary role to observational constraints 
arising from, for example, performing a PPN expansion (see for example |28| , |2!j| ). 

In Section 0, we begin an examination of the model on a flat spacetime, and show 
that the constraint imposed by A re-introduces the spatially longitudinal mode of 
the vector field as a physical field, giving it a propagation speed that depends 
nonlinearly on the vector field itself. We explicitly determine this velocity, and 
show that for spacelike ether models it is infinite near configurations with Aq = 0. 
In Section || we show that the linearization of these models does not correctly 
reproduce this result, and in Section we demonstrate that the field equations are 
also singular at this point. Also in this section, we show that the the field equations 
for timelike or spacelike ether models, will tend to form shocks. 

We also show, in Section [|, that these models generically do not have a classical 
energy that is bounded from below regardless of the signature of the ether, and 
in Section |^ we briefly consider the implications of adding additional terms to the 
action, and argue that such terms are not likely to help matters. Rather than 
proceed to a covariant or Hamiltonian analysis of the model in a curved spacetime, 
in Section ^| we consider a plane-symmetric spacetime, and reproduce the same 
propagation velocities as were observed in a flat spacetime. We conclude in Section^ 
with some overall comments. 

2. Flat Spacetime 

We begin by considering the model (|l]) in flat spacetime, since it is vastly simpler 
to illustrate some important dynamical aspects of the model without dealing with 
gravitational interactions. That we can do so and obtain results that are applicable 
to curved spacetimes is not a trivial issue though, as we shall describe more fully in 
Section |[ In any case, this sort of flat spacetime approximation (that is, so that one 
can transform g^ ss rj^ in a region of spacetime and consider small perturbations 
about this background) should be reasonable if g^ is in some sense the physically 
meaningful metric of spacetime. 
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The Lagrange multiplier field A enforces the constraint: 

A 2 = A 2 - A 2 = I, (2) 

where the sign of I indicates whether the vector field A^ describes a timelike (I > 0), 
null (I = 0) or spacelike (I < 0) ether. Here we are writing the spacetime covector 
field in terms of a scalar and spatial vector as: A^ — (A Q ,A), and will use, for 
example A - B = J2iAiB t and d 2 = Variation of (Q) with respect to A^ 

leads to the (nonlinear) field equations: 

d u F» v - XA^ = (3a) 

Contraction of this with respect to A^ and using (§) yields 

IX = A^ V F^ - A^J", (3b) 

which determines A provided that I ^ 0. Taking the divergence of (|a|) yields the 
constraint 

d„[\A»] = -d^. (3c) 

The first interesting feature of these models is that they possess multiple light 
cones, or multiple characteristics, indicating that there are fields in the theory that 
propagate at different speeds. Since this is not obvious by cursory examination 
of ([3a]), we will explicitly determine these propagation velocities. 

2.1. The Reduced Timelike and Null Field Equations. In both of these cases 
(I > and I = respectively) we can use (Q) to determine A$ since the only way 
that Aq can vanish is if A u vanishes identically. 



In the timclike case (3_b) will determine A, and the three remaining equations 
in (|3a| ) should determine A. To examine these equations we introduce a vector 
s M such that s^A^ = 0, which is a four vector = (so, s) with arbitrary spatial 
components Si, and timelike component determined from 

s = ft ■ s, (4) 

where we use the definition 

a=j- (5) 

^o 

Summing ( pa] ) into s M and using (^) to write what results in terms of Sj gives 

s^d v F^ = Si [ ni d a F 0a + 3 a F la ] = Sl (n l J° + J 1 ). (6) 

Since Si is arbitrary, the vector in square brackets is the three field equations that 
we are seeking: 

mdjFoj + d t F m - djFij = -{n l J° + T). (7) 

The same equation is arrived at in the null case by solving the zeroth component 
of (||a|) for A (this is a true algebraic equation for A since it involves no second-order 
time derivatives of any field) and substituting it into the spacelike equations. 
Re-written in terms of the vector potential, ([?]) becomes 

d 2 A t - dAAa - nAdjA^ - B 2 A, + n t d 2 A a + ii,i) r \, = n t J° + J\ (8) 
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These equations can be written solely in terms of A using the derivative relations 
that follow from 

d t A = ft ■ d t A, d t A Q = ft ■ d t A, (9) 
which leads to a second-order, quasi-linear partial differential equation of the form: 



where 



Lij(A k , dl d t d k , dkdfiAj + bi(A k ,d t A k , d t A k ) = 0, (10a) 

Lij = [Sijdt - iijdtdi - riidtdj - S^d 2 + n i n j d' 2 + didj] , (10b) 
b t = A^ 1 n i (Ski ~ n k ni)djA k djAi - A~ 1 (S jk - njn k )diAjd t A k . (10c) 



While the lower-order derivative terms (10c) are important to determine the partic- 



ular form of the solutions of (|10D, in this case they are irrelevant when determining 
the characteristics (see, for example, the discussion in [[30|). For this second-order, 
quasilinear system, the characteristic polynomial (see in particular the discussion 
in Section VI. 6 in |3l| ) may be found by making the replacement di —> ipi and 
dt —* ipa in L^, and computing the determinant 

det(L tf ) = (p 2 - f) 2 (p -n-p) 2 =0. (11) 

(This is most easily accomplished performing a rotation so that n = (m, 0, 0) which 
we can do since we are working pointwise-the results can then be re-assembled to 
give the general result.) This shows us that the initial value problem is well- 
posed provided we are imposing initial data on non-characteristic surfaces, where 
is this case the condition for a hypersurface to be non-characteristic depends on 
the vector field Ai. This is not particularly unusual, as exactly this situation 
also occurs in general relativity where the condition that data be imposed on a 
spacelike hypersurface involves the metric tensor. While ( |To| ) is not diagonal, and 
so the discussions in, for example |3^| , do not apply directly, it may be made 
so following the procedure given in |)3lf . 

We see therefore that there are characteristics defined by 

p = ±\p\, and p =fi-p, (12a) 

two of which are the usual null propagation conditions, and the third depends 
on nonlinear contributions from A. We therefore see quite clearly why this mode 
was not so well understood when linearizing about a flat background with A = 
p6| : in the linearized limit we would have a po = 'characteristic', which could 
be interpreted as a non-propagating field (that this mode should propagate was 
recognized in Jl6|). 



Before considering the spacelike case, we note that we can write (12a) in the 
covariant form: 

p 2 = 0, and A^tf- = 0, (12b) 

and so it is not too surprising that we will find exactly the same result. As we will 
see in Section |], and later on in a curved spacetime in Section [z| the field equations 
associated with the new characteristic surfaces have a tendency to form shocks. 
Note that although we have performed this analysis in three spatial dimensions, 
the result should hold in any number of spacetime dimensions. 
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2.2. The Reduced Spacelike Field Equations. It requires a bit more effort to 
show that dl2] ) also holds in the spacelike case. To begin, we write = (Aq, Aa) 
where o is a unit vector in the direction of A and A is the magnitude of the 

spatial component (so A = and a = A/ A). The normalization condition is 

Aq — A 2 = I < 0, which can be solved for A (this requires that A > 0, and it is now 
possible that Aq = 0). 

We again introduce a vector satisfying s^A^ = 0, but in this case we have 
so A) = As ■ a, and therefore the vector s must have the form 

s = Sao, + st, where st • a = 0. (13) 

Again summing s M into ( |3a|) gives only three equations, one proportional to sq and 
the other two proportional to St' 

s (d v F 0v + ^■a i d u F b '} + s Tl d„F l » = s (j° + ^J*) + s TiJ\ (14) 
which can be written as 

dvF Q V + ^L a . dvF iv = J° + ^-a 3 J\ (15a) 
(<% - a,u ,)i)J i" = (6 id - n,n,:.r. (15b) 
Multiplying the second by Aq/A and adding the first multiplied by a*, we have 

OidyF^ + ^d v F iv = a t J° + ^ J\ (16) 



which are three equations that correctly reproduce (|15|). 

Since ( |l6| ) is equivalent to (0) up to a multiplicative factor Aq/A, the analysis 
that we performed to find the characteristic surfaces will proceed as before, and 
lead to (|l^). The analysis has properly taken care of the fact that here A cannot 
vanish and Aq can, and we see from ( |j"2| ) that in the Aq — > limit the speed of 
propagation of one of the modes goes to infinity. 

3. Perturbative Analysis 

We consider a generic background A^ and perform a linearized expansion: A^ — > 
A^ = SA^, A = A + <5 A and J M — > J M + SJ^. From the expansion of (|^) we see that 

A^SA^ = 0, (17a) 

regardless of the sign of I. Then provided that A ^ we can proceed in a similar 
manner as in Section |^, which leads to the perturbation equations 

djdFij - dtSFio - UidjSFoj + A(% - n.n^SAj = 5 J 1 + ni SJ°. (17b) 

From this we see that this has a principle part that is equivalent to (Q) , and therefore 
perturbations will have characteristics determined by ([l2|). This result is redundant 
given that we had already determined the characteristics in the previous section, 
but may be looked upon as an alternative derivation. 

That ( |l7b| ) is correct unless Aq = 0, should be a signal that there is something 
special about such configurations. What we are going to show is that performing 
a linearization about such background solutions does not produce linearized equa- 
tions that approximate exact solutions of the field equations. This is a Cauchy 
instability [ p3[ in the exact theory, which shows up as a linearization instability in 
the perturbative calculation (see, for example pll |35|). In this section we will show 
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that the linearized equations only have p 2 — characteristics, contradicting the 
exact results (|l2|). We will also show in Sections |], || and]?] that the field equations 
are singular in the Aq — > limit. 

We must first consider background solutions with Aq = 0. Note that in order to 
write the field equations in the form ( fl6| ) we multiplied (15b) by Aq, and therefore 
taking Aq = in ( |l6| ) will not result in a full set of field equations. From ^) 
we must have A 2 = —I > 0, so we see that we only have two degrees of freedom 
remaining. The zeroth component of ( |3a| ) results in 

dtdiA = -J°, (18) 

and (|3l]) determines A as: 

IX = AiididjAj - B 2 A l + d 2 Ai) - AiJ\ (19) 

Since this determines A we can use this to write the remaining two spatial compo- 
nents of ( |3a"| ) as 

Pij {d 2 A j + d 3 d k A k - d 2 A, ) = P id F , (20a) 

where the projection tensor is defined as 

p ij = S v +l~ 1 A l A j . (20b) 

Note that we have three equations for the two remaining functions in Aj. This is 
not surprising since we have imposed Aq = by hand, and we expect that this will 
only be possible with some restriction on the source J M . Rather than determining 
the general form of these restrictions, we will choose a particular solution with: 
J° = and J 1 = 0, which implies that A = and allows a constant Aj. 

Performing the linearized expansion about this solution, from (17a) we find that 
AiSAi = 0. The timelike component of the linearization of ( |3a| ) leads to 

dfSAo - d t di5Ai = 6J°, (21a) 



and summing Ai into the linearization of the spacclike components of (17a) yields 
an equation for SX: 

ISX = -Ai(d t 6Fi - djSFij) - Ai6J\ (21b) 

Using this to remove SX from the remaining perturbation equations results in 

Pij {d t 8F l0 - djSFij) = -PijSP . (21c) 

Assuming that Aj is a constant vector pointing along the x-direction: Aj = 
(a, 0,0) where a 2 — —I, then (17a) implies that 6A X = 0. It is then a fairly 
straightforward matter to re- write the field equations that result from ( pl| ) as: 

(d 2 - d 2 )dlSA = {d 2 - d%)6J° + d t (d y SJ y + d z SJ z ), (22a) 

(d 2 - &)dl(d y 8A v + d z 8A z ) = d 2 x {d y Sjy + d z SJ z ) + d t (d 2 + d 2 )SJ°, (22b) 

(d 2 - d 2 )(d y 8A z - d z 5A y ) = d z 8J v - d y SJ z . (22c) 

This system contains three modes that all propagate at the speed of light, with no 
sign of the mode that propagates along the p^A^ — characteristic. 

These linearized equations would therefore lead one to believe that small fluctu- 
ations about a field configuration with Aq = all propagate at the speed of light. 
We know from the exact field equations that this is most definitely not the case, 
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and that there is a degree of freedom that has an infinitely large propagation ve- 
locity in the limit that Aq = 0. It should be reasonably clear that this is happening 
because we are attempting to perform a perturbation analysis about a degenerate 
field configuration-in fact our initial data surface is a characteristic surface for the 
exact field equations. 

This suggests us that the expansion performed in Q is suspect. It is, of course, 
possible that other models will result in the linearized equations that they use, but 
it is clear from what we have seen here that the linearized equations that they 
derive do not approximate solutions to the exact field equations. Lest one think 
that this sort of thing would not happen when the model is coupled to gravity, in 
Section |?j we show that the same characteristics (|l^) appear in a curved spacetime. 

We also see that this does not happen in the null and timelike cases, but this 
does not mean that there are not other issues to deal with. We will show in the 
next section that shock formation should be expected, and in Section 0, that the 
(classical) Hamiltonian of these models will not be bounded from below. 

4. Shock Formation 

Now we wish to examine some properties of nontrivial solutions that will display 
propagation along the p^A^ = characteristic. One such example is that of a plane- 
symmetric vector potential: A^ = (Ao(t, x), At(t,x),0, 0) (and correspondingly for 
the source), since these solutions are expected to propagate along the x-direction, 
along which A also points. What we will discover (for a non-null ether) is that the 
field equations reduce to a form that can be written as a conservation law and is 
generically expected to form shocks. 

From ([5a]) we find the two field equations: 

d 2 x A Q - d t d x At - XA = J°, (23a) 

d%Ai - d t d x Ao + XAi = J 1 , (23b) 

and taking a linear combination that removes A we find: 

A l {d 2 x A - d t d x At) + A (d 2 t A 1 - d t d x A Q ) = A X J° + A j\ (24) 

Using (H), this can be written as an equation determining either Aq or At, both 
forms of which are: 

(d t - ^d x ) (d t At - ^d x At) = + J\ (25a) 

(^dt - d x ) (^-d t A - d x A Q ) = J° + ^ J 1 . (25b) 

In either case we see that in a region external to the sources we need only consider 
a homogeneous, first-order, partial differential equation, solutions to which are also 
solutions when the source is non- vanishing. 

For a null ether we have Aq — ±At, and ( |25a| ) leads to 

d t Ai T d x At = 0, (26a) 

with solutions: At = Ai(t±x). Throughout we will assume that Aq > (the ether 
is forward-pointing in time), and so disturbances in the ether field propagate in 
the opposite spatial direction to which the ether field itself points. This last fact 
is a generic feature, consistent with the characteristic surface defined by p^A 1 - 1 = 
in (0). 
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For a timelike either we can solve (|J) for Aq as: Aq = ±y/ A\ + 1, and so we 
consider solutions of: 



8 t A lT 7 =^= d x A 1 = a t A 1 Td x J^+l = 0, (26b) 
V A 1 + I 

where the second form is written as a conservation law. This type of equation is 
considered in many standard texts on fluid or gas dynamics (see for example [ (30[ 
[36|), from which it is evident that, for example, a Gaussian fluctuation of A\ about 
a constant background Aq will inevitably lead to shock formation. The value of A\ 
is constant along the characteristics, which are straight lines determined from the 
initial value of A\ by: x(t) — k(xo)t + Xo where k(xo) — ^A\{xq) / ^ A\{xq) + I and 
A\{xq) = Ai(0, xq). Note though that since I > we have \k\ < 1, and disturbances 
propagate at a velocity slower than that of light. 

For a spacelike either we have a similar situation. Solving (Q) for A\ and us- 
ing ( [25b| ) to determine Aq, we find (this can also be written as a conservation 
law): 



dt A T - d x Ao = 0, (26c) 

which shows that disturbances not only generically form shocks but that they also 
travel faster than light. We also see quite clearly that the limit Aq — > does 
not lead to a smooth linear limit. Writing the characteristics as above, we now 
have k(xo) = ^yAg — 1/Aq, and regions with very small A have characteristics 
that approach infinite velocity. Since disturbances in Aq cannot be assumed to be 
exactly constant over an extremely large region, shock formation is to be expected. 

Although we have explored the system in a plane-symmetric situation (which we 
shall explore further in Section Q in a curved spacetime), we have found that an 
identical result holds for spherically-symmetric solutions. Although a spherically- 
symmetric ether is difficult to motivate, it does show that the existence of shocks 
is not an artifact of the specific solution that we are considering. 

At this point it is worthwhile saying a few words about a restricted gauge invari- 
ance that the system possesses when the external source coupling is conserved (that 
is, when d^J^ =0). In this case the action (|l|) is invariant under the gauge trans- 
formation A^ — > A^ + d^9 provided that 9 satisfies: g^ v {A fjL + d^9)(A u + d u 9) = I. 
In the present model, this simplifies to (solving for d t 9): 



d t e = -A ±^Al + d x e(dJ + 2A 1 ). (27) 

It can be shown that all of ( p6| ) are invariant under this transformation, so that A^ 
may be transformed to zero by an appropriate choice of gauge. 

This should not be too surprising, since in this case the model is identical to a 
Maxwell field in the nonlinear gauge: A 2 = I. Therefore, if the ether field is only 
coupled to a conserved current then the model is gauge-equivalent to a Maxwell 
field, and there would be little point working in the nonlinear gauge (||) when 
some more well-behaved linear gauge (such as the Lorenz gauge d^A^ = 0) would 
equivalent and far be simpler to work with. The intended interpretation of the 
model is that A^ should represent an ether field, and therefore it will couple to 
matter fields as a timelike metric representing the preferred frame: g^ u ~ A^A V [ fl6| 
p7[ , and so coupling to a conserved current is neither expected nor desired. For the 
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ether to have observable consequences it must couple to a non-conserved current, in 
which case the p^A^ = characteristic is physical, and shock formation is expected. 



5. Canonical Analysis 

The existence of multiple light cones and shock formation in a theory of the 
ether is perhaps not sufficient reason to discard the model. The theories that we 
are considering here have another dangerous feature though: they do not have a 
classical Hamiltonian that is bounded from below in flat spacetime. This analysis 
should be enough to convince the reader that in any (stable) region of spacetime 
that can be approximated by a flat spacetime, these theories should be generically 
unstable once quantum or other dissipative effects are included. 

In a quantum scenario, the standard argument is that quantum fluctuations can 
result in the creation of real (as opposed to virtual), negative energy particles, since 
no energy need be 'borrowed' from the vacuum to create them. Thus the quantum 
system will cascade into states of increasingly more negative energy. While it is 
possible to claim that these negative energy states are not realized in a classical 
model (if the initial data has finite energy it will remain with the same finite energy) , 
vacuum fluctuations in the quantum theory seem unavoidable. 

The analysis is fairly standard, so we will be brief. Expanding the Lagrangian 
from (fll) we find 



L = d 6 x 



^(dtAt - dtAo) 2 - h,F l3 ) 2 + ^\(A 2 -A 2 -l) + A J° + A-J 



(28) 



and we see that we can treat Aq and A as Lagrange multipliers, and the conjugate 
momentum for A is 



C T 

P * = "FTT" = 9 t A i ~ d * A o- (29) 



The Hamiltonian is then 
H = I d 3 xP i d t A l - L 



= I d 3 :, 



\P 2 + P- VA + l(Fij) a - l -\{Al -A 2 -I)- A J° -A- J 



(30) 



and we will assume that the fields fall off fast enough as r — > oo that no surface 
contributions are required. 

The condition imposed by A is (|j), which can be used to determine the Lagrange 
multiplier Aq as 



A ^±\]A* + l. (31) 

This is a delicate point since, as noted before, if I < there is no guarantee that 
Aq 7^ 0, but since (||) involves the the Lagrange multiplier Aq it cannot be considered 
as a constraint. The condition imposed by Aq is 

^L^-V-P~XAq- J°«0, (32) 
oAq 
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so that since Aq is already determined, this determines the Lagrange multiplier A 



V • P 4- 7° 

A«T , F + J • (33) 

Va 2 + i 



as: 



Evaluating the Hamiltonian (|30|) "on-shell" (replacing the Lagrange multipliers 
by their phase space representation using ( |3l| ) and (|33|)) in a region where the 
source terms vanish, we find the remaining contributions: 

■^i^M + k).], (34) 



Hea d 



j3 . 



where we have used c^Ao = ±(AjdiAj)/y/ A 2 + 1. To show that this can be ar- 
bitrarily negative on phase space, consider configurations derivable from a static 
scalar potential as A* = di(f>(x). The conjugate momentum to A is Pi = —diAo, 
and evaluating the Hamiltonian density yields: — (1/2) (diAo) 2 , which can be made 
arbitrarily large and negative by choosing <f> to be a Gaussian with large amplitude. 
Note how this has occurred: the imposition of (0) has caused (|32|), which would 
have been the Gauss constraint in Maxwell's theory and led to a positive-definite 
term in the Hamiltonian, to now become an equation that determines a Lagrange 
multiplier, and the Hamiltonian is indefinite. 

We can also consider Hamilton's equations, which appear fairly innocent when 
written in terms of the Lagrange multipliers: 

At = {A, H} = ^-=Pi + d t A « P ± ' (35a) 
sp i v A 2 + I 

SH 

Pi = {P h H} = -— = -djFy - XAi + r 



dj F i3 ± Ai — / = +J\ (35b) 



V • P + J' 

Va 2 + i 

where in the second form of these ( |3l| ) and ([53|) have been imposed. Note that as 
far as the initial value formulation is concerned the sign of I is largely irrelevant, 
but manifests i tself by allowing the possibility that Aq = in the spacelike case, 
in which limit ( |35b| ) is singular. This is another reflection of the fact that in this 
model Aq — is a special point, and we cannot expect that it may be reached 
smoothly (see also [p7|). 

6. Related Models 

Before demonstrating that all of the aforementioned issues also arise in the curved 
spacetime form of the model, we make a few comments on flat spacetime general- 
ization of (Q). 

If we were to introduce non-Maxwell kinetic contributions such as (V M A M ) 2 
into (jl|) as suggested in jf7|, we would perhaps have different characteristics, but 
all of the difficulties that we have been discussing here would be present. Since the 
model will have time derivatives of Aq in the action, Aq is no longer a Lagrange 
multiplier and the P ■ 8Aq term remains in the Hamiltonian, along with new con- 
tributions Pq j (2d) + PqB ■ A involving Po (the momentum conjugate to Aq). In this 
case since both Po and Aq are phase space coordinates, and it is a trivial matter 
to find field configurations on which the Hamiltonian is as large and negative as 
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one likes. In order for this theory to be equivalent to a Maxwell field written in a 
nonlinear gauge one would need a restricted gauge invariance that simultaneously 
left the (V ^A^) 2 term and the condition invariant, which seems unlikely. That 
this type of model is 'dangerous' in this sense has already been discussed in the 
literature p^j. The same reference also shows that the model considered as an 
example of a preferred- frame theory of gravity in |28| , has stability issues that go 
beyond the PPN analysis performed there. 

Similar comments apply to the possible addition of curvature-coupling terms such 
as R^ V A^A V . Since this would constitute an explicit derivative-coupling, altered 
characteristics are to be expected p2| ], and were observed in ]37| . In the latter 
reference the characteristic surfaces of the exact theory are not determined, and 
so it is not known whether the perturbation that is performed is degenerate in the 
sense discussed in Section ||. 

If we alter the term that imposes the constraint (^) to X(A 2 — I) 2 as suggested 
in then A will drop out of the field equations altogether, leaving (||) and 

d v F^ v = J M . Consistency of the field equations requires that the vector source is 
conserved: d^J^ — 0, and since the divergence equation is trivial there are really 
only four field equations. Since the coupling is therefore required to be gauge- 
invariant, this is a Maxwell-theory written in a nonlinear gauge defined by (^|). 
Therefore the n = 2 model of jl8) is just Einstein-Maxwell theory in disguise, with 
the Lorentz symmetry breaking as the result of a nonlinear gauge choice for the 
vector field. It is interesting that adding a self-coupling potential term of the form 
V(A 2 ) to this theory could have disastrous consequences, since even if the coupling 
current is conserved the divergence of the field equations lead to d^A^ = 0, which is 
a nontrivial equation determining a degree of freedom in A^, therefore there are five 
field equations for the four unknown functions A^. In contrast, introducing a self- 
interaction term to (|l]) merely induces the shift: A — > A + V'(l). This demonstrates 
that the inclusion (or not) of a mass term in a vector field theory can have rather 
significant consequences, contrary to the statement made in |37| , Section 5.4]. 

We should mention that the vector field-based model that we have introduced [[n] 
is of a different class than those discussed here. Because we are imposing no addi- 
tional constraints in addition to those already present in an Einstein-Proca model, 
we do not expect the issues discussed in this paper to arise. On the other hand, 
since the vector field is not guaranteed to be timelike in these models (although 
this turns out to be the case in a cosmological setting), it is not so obvious that we 
can interpret the vector field as representing a dynamical ether fluid in general. 



7. Plane-Symmetric Spacetimes 

We will now examine the model in a particular curved spacetime, demonstrating 
that the characteristics (|l2|), existence of shocks, and problems performing a linear 
expansion near Aq = all survive. This is a curved spacetime generalization of 
the plane-symmetric solutions considered in Section |?], which is general enough 
to contain the static background solutions that were considered in |Q. We can 
therefore examine some nonlinear solutions 'near' this spacetime. 

Working now in a curved spacetime, from (|l|) we find the field equations for A^: 



grA„A v = I, V„F"" - \A» = 0, 



(36a) 



12 



M. A. CLAYTON 



and the Einstein equations: 

<V = -\(F m F v a - \g^F 2 ) + X -\A^A V . (36b) 
The spacetimes we are interested in have an 'ether wind' along one direction: 
= (A o (t,x),A 1 {t,x),0,0,... ,0), (37) 

consistent with translation and rotation killing vectors in a flat n-dimensional hy- 
persurface, with plane-symmetric metric: 



. _ fgab(t,x) 

\ -a i ' n {t,x)5 lJ 



(38) 



where a, b, ... € {1, 2} and . . . € {3,4, . . . ,n+ 2}, so that n is the dimension of 
the flat subspace. We use coordinates {y 1 } in the fiat hypersurface and t and x in 
the 'curved' directions. We will also use the remaining coordinate invariance in the 
t — x plane to choose the conformally-fiat metric 

gab = e^Vab, (39) 

where rj a b is the two-dimensional Minkowski metric: r\ a \, = diag(l, —1). 

Reducing the action results in (dividing out the infinite volume J d n y in the flat 
plane) : 



Sred = I dtdx 



-2an ab d a ad b (f> - — ^-r, ab d a ad b a 
n 



\a 2 e^F^ + \\a 2 (jf b A a A b - Ie*)l , (40) 



where 

F 01 = d t A x - d x A . 
From this we find the field equations 

{dt- d l)c? + c?(l e -*F* 1 - 1 -lM 

(^_ a ^ + 4(n-2)i 2 _ dl)a-\e 
n a 2 
Al - A\ - le* = 0, 

d x [a 2 e-*F 01 ] -2Xa 2 A Q = ^ 

a t [a 2 e-^F i] - 2\a 2 A l = 0, 

which are identical to what we would have found directly from the field equa- 
tions ([56]). The linear combination of the last two of these from which A cancels 
gives 

A d t [a 2 e-*F 01 ] - A x d x [a 2 e^F 01 ] = 0, (43) 

and we see that as expected, the p^A^ — characteristic has survived. 

Although ( (43| ) results in a 2 e~*i<oi = constant, along these characteristics, it 
is difficult to solve the remaining equations. We will instead consider the simpler 





(41) 


= o, 


(42a) 




(42b) 




(42c) 




(42d) 




(42e) 
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class of solutions with Fqi = 0, which also implies that A = 0. In this case we find 
(writing the equation in both forms analogous to (|25|)) 

MdtAx - A 1 d x A 1 = \le*d x cf>, (44a) 

A Q d t A - A 1 d x A = \le*d t 4>. (44b) 

Both of a and <fi will be functions of either t + x or t — x, and for the the solutions 
4> that are either right- or left- moving {(j) — (f>(t =p x)), we find (again assuming 
for forward-pointing Aq > 0) that A\ propagates in the same direction as (f>. The 
equations (Q) then become (<p' — dt4>{t =F x)): 



dtAo T d x A = lie*?-, (45b) 

At) An 



lo - m 

and the gravitational wave 4* is unaffected by the presence of the ether. 

We have the following scenarios for the different types of ether: For a null ether 
all fields are decoupled, and all travel at the speed of light. For a timelike ether, 
the gravitational wave 4> acts as a source in (45a), disturbing the ether. In this 
case disturbances in the ether travel slower than 4>, and we generically expect shock 
formation in the ether. For a spacelike ether we have a similar situation, but (45b) 
tells us that the passing of the gravitational wave produces ether fluctuations that 
travel faster than light, and therefore faster than the gravitational wave disturbance 
itself. Similar to the flat spacetime case discussed in Section ||, the limit Aq — ► 
does not allow a linearization that represents the behavior of A Q as determined by 
the exact field equation (45b). 

It would be interesting to explore more general solutions to ( f42| ) to find out how 
the back-reaction affects the propagation of gravitational waves. Although for the 
simple _Foi = solutions that we have considered it is a fairly simple matter to see 
that it is the passing of the gravitational wave that is disturbing the ether, in the 
more general case this is not obvious. Since this would probably require numerical 
analysis and since the speed of light becomes very large as Aq — > 0, satisfying the 
Courant-Friedrichs-Lewy condition could be a delicate business. 



8. Discussion and Conclusions 

Although the nonlinear constraint imposed in these models is algebraic (A 2 = I 
for some fixed I), we have shown that it nonetheless alters the causal evolution 
of the vector field. This is well-known when working in non-covariant gauges in 
Maxwell electrodynamics, where the resulting light cones should have no physical 
consequences. In the case at hand, the (co-)vector field A^ is intended to repre- 
sent a fundamental ether frame, giving a dynamically-determined, timelike, null or 
spacelike direction in spacetime. If this is the case, then it should not be obser- 
vationally equivalent to a Maxwell field, the above constraint will have nontrivial 
consequences, and therefore the multiple light cones that we have shown exist in 
such models will be observable. 

We have also shown that these models have a classical Hamiltonian that is not 
bounded from below. Although not necessarily a concern in classical conservative 
systems, we have to expect that any proposed quantization should seriously address 
the existence of a stable quantum vacuum. We have argued that things become 
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worse in this respect when terms like (V^^ 1 ) 2 are introduced into the action, since 
then all four components of would have conjugate momenta, and we have even 
more freedom to find classical field configurations of arbitrarily large, negative en- 
ergy. In addition, given the analysis of propagation speeds in linearized analyses of 
similar models [ j37| , a simple light cone structure is not to be expected. Introducing 
curvature-coupling terms into the action cannot be expected to alleviate matters 
since the derivative-coupling that results is known to introduce altered propagation 
speeds for different fields p^ ], which has been explored in very similar vector field 
models in p7| . 

One of the main messages of this work is that it is important to realize that a 
linearized analysis does not necessarily tell the whole story (i.e., a theory may have 
a 'linearization instability'). In the case of a spacelike ether we saw that it is quite 
possible to formally expand about a degenerate configuration, and although we 
know from a more complete analysis that there is a field with an infinite propaga- 
tion velocity, this never appears in the linearized analysis. The canonical analysis 
has also shown that the Hamilton equations became singular at these solutions, 
and we therefore cannot expect that the resulting linearized equations to bear any 
resemblance to solutions of the exact field equations. This is precisely what was 
noticed in earlier work on vector fields coupled to gravity |27j , and shows that the 
analysis in jl8| is suspect. The same concerns may arise when considering the lin- 
earized analysis of vector models in | f37| , since the models are derivatively-coupled 
and the characteristics have not been determined. 

In any case, the quantization of any of the models that we have discussed here is 
delicate, since the classical theory does not possess simple characteristics, fields in 
model are shock-forming, the classical energy is not bounded from below, and the 
model possesses linearization instabilities. 
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